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Abstract 

Active liquid crystals or active gels are soft materials which can be physically realised e.g. by preparing 
a solution of cytoskeletal filaments interacting with molecular motors. We study the hydrodynamics of an 
active hquid crystal in a slab-like geometry with various boundary conditions, by solving numerically its 
equations of motion via lattice Boltzmann simulations. In all cases we find that active liquid crystals can 
sustain spontaneous flow in steady state contrarily to their passive counterparts, and in agreement with re- 
cent theoretical predictions. We further find that conflicting anchoring conditions at the boundaries lead to 
spontaneous flow for any value of the 'activity' parameter, while with unfrustrated anchoring at all bound- 
aries spontaneous flow only occurs when the activity exceeds a critical threshold. We finally discuss the 
dynamic pathway leading to steady state in a few selected cases. 
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INTRODUCTION 

Active viscoelastic gels and active liquid crystals are a new kind of soft materials iH S, 0, 
[sj]. They differ from their passive counterparts in the crucial respect that they continuously burn 
energy, e.g. from chemical reactions, which drives them out of thermodynamic equilibrium even in 
steady state. Examples of passive viscoelastic fluids are polymeric fluids, rubber, molten plastics 
etc. Examples of active gels are cytoskeletal slabs biomimetic systems iMor solutions of 

molecular motors such as dyneins or myosins and microtubules or actin fibers [ftTjlsiR 10. 11]. 
In these solutions the molecular motors interact with actin or microtubules by exerting forces or 
torques on them, creating motile and transient cross-links etc. It is the presence of the molecular 
motors, which are able to use chemical energy to perform mechanical work, which renders these 
viscoelastic polymeric fluids 'active'. The source of chemical energy in such active gels typically 
comes from ATP hydrolysis (there is estimated to be ~ 15 ksT worth of biochemical energy 
from each ATP hydrolysis reaction). It has also been suggested that the equations of motion of a 
continuum active gel may explain some macroscopic aspects of motility in bacterial colonies, such 
as swarming (i.e. the collective behaviour of a large number of bacterial swimmers), elastotaxis 
etc. iBB]. 

In general, the modelling of active soft matter presents novel fundamental challenges to theo- 
reticians, due to the non-existence of a thermodynamic free energy. Furthermore, the elastic and 
rheological properties of these materials may be anticipated to be strongly affected by their non- 
equilibrium nature Tactivity'). E.g., we now know that spontaneous flow can exist in steady state 
in active materials [llil2|], in sharp contrast to what happens in their passive counterparts, where any 
non-zero current is disallowed in equilibrium by momentum conservation and the Navier-Stokes 
equations. Furthermore, many biological gels, such as actin and neurofilament networks, thicken 



when sheared 111211 . which is the opposite of the typical behaviour of viscous polymeric fluids, 
molten plastics etc., which flow more easily as the external forcing driving their flow increases. 

Our goal here is to study numerically the equations of motion of an active liquid crystal sample 
confined between two parallel plates in a quasi- ID geometry, focussing on the effects of different 
anchoring conditions for the director field (the active analog of the director field), and of the fluid 
behaviour deep in the active phase. Numerical methods are needed in this, as the insight which 
is possible to obtain analytically is limited although important. For instance, it has been possible 
to theoretically predict M\ by analytic methods the existence of an activity-induced transition 
between a passive phase with no local flow to an 'active' phase, in which the fluid is flowing at 
steady state. However, as the techniques used thus far essentially allow a stability analysis of 
the problem in selected geometries, it is not possible to predict e.g. the flowing and director field 
patterns deep in the active phase, away from the transition to spontaneous flow. Furthermore, using 
a numerical method allows us to consider more detailed forms for the equations of motion. Thus 
here we extend the vector polarisation model considered analytically in Ref. [HiQ], but a tensorial 
equivalent (similar to the one introduced in Ref.JsiH]), which allows us to study situations in which 
the order parameter is not constant throughout the sample (as occurs in passive liquid crystals e.g. 
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in shear banding and close to a defect core). 

We will show that the existence of elastic frustration in the sample rising from the anchoring at 
the boundaries is important in determining the steady state solution, just as it is in selecting director 
field patterns of a passive liquid crystal under an imposed shear. Here we focus on two examples 
of boundary conditions. With normal anchoring, we find a transition between passive and active 
phase at a (size-dependent) non-zero value of the 'activity' parameter. With conflicting anchoring 
(we choose the anchoring leading to the hybrid aligned nematic - HAN - cell in conventional 
passive liquid crystals), on the other hand, the behaviour is quite different and the threshold is zero, 
i.e. any active liquid crystal in a HAN geometry spontaneously flows in steady state. Increasing the 
activity (going deeper in the active phase), results in other 'transitions'. In the normal anchoring 
case we observe banded flow, with the number of bands increasing as we go deeper into the active 
phase, while in the HAN cell geometry a Poiseuille-like flow sets up, and changes direction upon 
increasing the activity, resulting in vastly different elastic deformations of the sample in steady 
state. We also follow the dynamic pathway leading to steady state in selected cases and find that 
transiently the director field may be transiently asymmetric although it is eventually symmetric 
in steady state. Furthermore we observed some dependence of the steady state we find in our 
simulations on the initial conditions, hence on the history of the system. 

This work is structured as follows. In Section II we introduce the tensorial model we consider 
and the equations of motion we want to solve. We also briefly show how the lattice Boltzmann 
algorithm of Ref. [|l3ll may be generalised to solve these equations of motion. In Section III 
we report our results for the flow and director field patterns in steady state in the case of normal 
(Section III A) and conflicting (Section III B) boundary conditions at the sample walls. In both 
cases we assume translational invariance parallel to the walls so that the numerical problem is 
effectively one dimensional. We also discuss the dynamics leading to steady state in a selected 
case with normal anchoring (Section III C). Finally, Section IV contains our conclusions. 



MODEL AND METHODS 

We introduce a Landau-de Gennes free energy to describe the equilibrium physics of the active 
liquid crystal. This is a function of a tensorial order parameter, Qap, whose largest eigenvalue 
gives the strength of local order. The free energy JF is a sum of two terms, J^i + J^2- The first one 
is a bulk contribution. 



•^1 - -^(1 - ^)<5q/3 ^QafiQp-fQya + -^{Qlpf, (1) 



while the second, J^2 = K (d^Qais)'^ /2 is a distortion term lIlSll . Hereafter Greek indices denote 



Cartesian components and summation over repeated indices is implied. K is an elastic constant. 
We have taken the one elastic costant approximation and neglected the spontaneous splay co- 
efficient, which is equivalent to having infinitely strong anchoring at the boundaries. Adding a 
spontaneous splay does not change the physics qualitatively as argued in In the free energy 
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equation, Aq is a constant and 7 controls the magnitude of order. The anchoring of the director 
field on the boundary surfaces is along x or x, see Fig. 1. 



The equation of motion for Q is then 1II6I1 



{dt + u- V)Q - S(W, Q) = TH + AA/iQ (2) 

where F is a collective rotational diffusion constant, A/i is the chemical potential gained via ATP 
hydrolysis 111], and A is an 'activity' parameter of the liquid crystal. 

Eq. (I2I) is valid to first order in Afi The first term on the left-hand side of Eq. ^ is the 
material derivative describing a quantity advected by a fluid with velocity u. This is generalized 
for rod-like molecules by a second term 

S(W,Q) = (eD + ^)(Q + I/3) (3) 
+ (Q + I/3)(eD-a;)-2e(Q + I/3)Tr(QW) 

where Tr denotes the tensorial trace, while D = (W + W^)/2 and u = (W — W^)/2 are the 
symmetric part and the anti-symmetric part respectively of the velocity gradient tensor Wap = 
dpUa- The constant ^ depends on the molecular details of a given liquid crystal. In our context 
increasing ^ tends to disfavour tumbling and favour aligning of the director in the material. The 
molecular field H in Eq. [His given by H = -g + (I/3)Trg. 

The fluid velocity, u, obeys the continuity equation and the Navier-Stokes equation, 

p{dt + Updp)Ua = disiUap) + ^dp{daUp + d^Ua (4) 

where p is the fluid density and r] is an isotropic viscosity, and liap = Yi^^^^^'^ + E^'^^'^*'. is 
the stress tensor necessary to describe ordinary liquid crystal hydrodynamics: 

-Qpassivc _ _ p^^^^ _[_ 2^{Qap + -6al3)Q-ytH^^ (5) 

Pq is a constant in the simulations reported here. The active term is given by = CAiiQai3, 

where C is another 'activity' parameter. This equation is also valid at first order in AyU 11}]. 

The model above reduces to the Beris-Edwards model for liquid crystal hydrodynamics for 
A = C = 0. For uniaxial active liquid crystals, the active liquid crystal director field n is well 
defined through Q^p = 3g {uaUp — 5q,/3/3) /2, where q is the largest eigenvalue of Qap. In this 
limit our model can be shown to reduce to the vectorial model considered in 

We solve Eqs. [2]and|4]by using a Lattice Boltzmann algorithm. Our algorithm is a generalisa- 



tion of the one introduced in Refs. [|l3l] and successfully used to study a number of problems in 
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the hydrodynamics and rheology of a passive liquid crystal. We now review it briefly, following 
the treatment in Ref. [13]. 

Lattice Boltzmann algorithms to solve the Navier-Stokes equations of a simple fluid, are defined 
in terms of a single set of partial distribution functions, the scalars fi{x), that sum on each lattice 
site X to give the density. For liquid crystal hydrodynamics, this must be supplemented by a second 
set, the symmetric traceless tensors Gi{x), that are related to the tensor order parameter Q. Each 



fi, Gi is associated with a lattice vector lll3h . We choose a 15-velocity model on the cubic 
lattice with lattice vectors: 

ef^ = (0,0,0) (6) 
ef^ = (±1,0,0),(0,±1,0),(0,0,±1) (7) 
ef^ = (±1,±1,±1). (8) 

The indices, i, are ordered so that i = corresponds to e^"-*, i = 1, ■ ■ ■ 6 correspond to the e^^-* set 
and i = 7,14 to the ef'^ set. 

Physical variables are defined as moments of the distribution function: 

P = X] /i, pUa = ^ fiCia, Q = ^ Gi. (9) 

i i i 

The distribution functions evolve in a time step At according to 

+ elAt, t + At)- fi{x, t) = ^ [Cf.ix, t, {/,}) + Cfiix + elAt, t + At, {/*})] (10) 

G,{x + e^At,t + At) - G^{x,t) = ^ [CGi{x,t, {G,}) + CG^{x + eiAt,t + At, {G*})] (11) 

This represents free streaming with velocity Cj followed by a collision step which allows the distri- 
bution to relax towards equilibrium. The C's are collision operator (detailed below) while /* and 
G* are first order approximations to fi{x + e-iAt, t + At) and Gj(x + At, t + At) respectively. 
They are obtained by using At Cf j(x, t, {fj}) on the right hand side of Eq. (flOl ) and a similar sub- 
stitution in Eq. (fTTI) . Discretizing in this way, which is similar to a predictor-corrector scheme, 
has the advantages that lattice viscosity terms are eliminated to second order and that the stability 
of the scheme is improved. 

The collision operators are taken to have the form of a single relaxation time Boltzmann equa- 
tion, together with a forcing term 

Cf,{x,t,{f,}) = --(/.(x,t) - /f(f,t, {/,})) +p,(f,t, {/,}), (12) 
^/ 

CG^{x, t, {G,}) = - — (G,(f , t) - Gt'ix, t, {G,})) + M,ix, t, {G,}). (13) 

The form of the equations of motion follow from the choice of the moments of the equilibrium 
distributions /f' and G^^ and the driving terms pi and Mj. The distributions f^"^ are constrained 
by 

fi'^ = ^' 5Z fi'^^^'^ = PUa, Yl f^^i^^iP = -H^™ + pUo,Up (14) 
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where the zeroth and first moments are chosen to impose conservation of mass and momentum. 
The second moment of f^'^ controls the symmetric part of the stress tensor, whereas the moments 
of Pi 

^Pi = 0, y^^PiCia = dpTaiS, "^PieiaCip = ^ (15) 

i i i 

impose the antisymmetric part of the stress tensor, which we have called Tap. For the equilibrium 
of the order parameter distribution we choose 

= Q, ^ Gfeia = Qua, ^ G-'^ei^ei/j = QuaUf^. (16) 

i i i 

This ensures that the order parameter is convected with the flow. Finally the evolution of the order 
parameter is most conveniently modeled by choosing 

J2 M. = rH(Q) + S(W, Q) = H, J2 = E ^')"" 

i i i 

which ensures that the fluid minimises its free energy at equilibrium. 

Conditions (fT4l)-(fT71) are satisfied, by writing the equilibrium distribution functions and forcing 
terms as polynomial expansions in the velocity. The coefficients in the expansion are (in general 



non-uniquely) determined by the requirements that these constraints are fulfilled (see Ref. [|l3h for 
details). 

The active contributions then alter the constraints on the second moment of the //s and the first 
moment of the G/s distribution functions as follows: 

^ M, = H + AA/iQ, ^ ft'ciac^p = -nP™ + puaup - (AfiQ^p. (18) 

i i 

Alternatively, we can input the derivative of the pressure tensor as a body force (as in Eq. [T5]) . This 
scheme rigorously ensures no spurious velocity in equilibrium [ isl, 1^ and has been employed in 
the following. 



RESULTS 

In this Section we report the results obtained with lattice Boltzmann simulations of the dynamic 
behaviour of a slab of active liquid crystal, sandwitched between two fixed plates at a mutual dis- 
tance L along the y axis (see Fig. [T]for axis definition). As a result our system is effectively one 
dimensional (the velocity and the active liquid crystal director fields depend on the velocity gradi- 
ent coordinate, y, only) We shall consider different kinds of boundary conditions, and demonstrate 
that these affect the hydrodynamics of the fluid considerably. 

The model parameters in the simulations reported here were typically fixed as follows: Aq = 
0.1, L = 100, K = 0.04, 7 = 3,^ = 0.7, T ~ 0.34, while the relaxation time of the lattice 
Boltzmann evolution of the density distribution function (see Ref. H for details) is 0.7. AU 
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preceding values are in simulation units. These can be straightforwardly related to a set of physical 
units by considering that one space and time simulation units correspond to 0.02 yum and 0.067 
fis respectively (for how to do this mapping see Ref. lilBII ). One can thus map our sample onto 
a physical slab of active liquid crystal of thickness 10 yum, while the parameters of the liquid 
crystalline fluids are K = 16.2 pN, 7 = 1 Poise (7 is the rotational viscosity lIlSll ). while the 
ratio between the Leslie viscosities as/a2 ~ 0.08. Note that a positive value of as/a2 means 
that the liquid crystal in the passive phase is flow aligning niSll : thus the parameter choice made 
here corresponds to a flow aligning liquid crystal in the passive (( = 0) phase. We have further 
fixed A = in our simulations - as this parameter has been shown to not affect the^xistence of a 
transition between a passive quiescent and an active phase with spontaneous flow [2]. The values 
of the remaining 'activity' parameter were on the other hand systematically varied. 

In what follows the director field (also called polarisation field in the context of active gels in 
Ref. fl) remains in the xy plane in the vast majority of cases. It is therefore useful to introduce 
an orientation angle 9, as the angle between the director field and the positive direction along the 
y axis. 

Note that, unless specified otherwise, we always initialised the system with zero velocity field 
and with a director field corresponding to the stable equilibrium one in the passive phase. How- 
ever, as this initial state is generally metastable, and as our lattice Boltzmann simulations do not 
include noise, we introduced a small perturbation - typically in the form of a small deviation in 
the orientation angle at y = L/2. This is enough to allow the system to escape from the metastable 
passive phase in which it is initialised. 



Active liquid crystals with normal anchoring 

We first consider an active liquid crystal with normal anchoring at its boundaries (Fig. la). 

For values of ( smaller than Q ~ 0.0015 the liquid crystal remains in the passive phase, whereas 
for C > Cc there is a finite flow along the x direction (Fig. [T]), i.e. the liquid crystal is in an active 
phase. Fig. [21 shows the director and velocity profiles first encountered in the active phase upon 
increasing the value of C- There is a net primary flow along the x direction, which exhibits a sharp 
maximum at the channel center, and is quite different from e.g. the common parabolic Poiseuille 
velocity profile which may be externally induced in a passive liquid crystal by exerting a pressure 
difference or a body force on it. The spontaneous flow symmetrically deforms the director profile 
and two 'bands' develop, in which the orientation angle 9 attains opposite values (close to that of 



the Leslie angle [11511 which is imposed by the local shear rate introduced by the spontaneous flow). 
These band flatten when the value of ( is increased - up to a second critical value, which for the 
parameters we use here is ~ 0.0015. 

A further increase in leads to a qualitatively different behaviour, as shown in Fig. [3l The 
velocity profile changes and now has a sinusoidal shape, so that it is zero at the centre of the 
channel and there is no net flow. For values of ( close to this second crossover the director field 
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exhibits an out-of-plane instability, much as conventional (yet flow tumbling) liquid crystals under 
strong shear. For large the director field is well aligned along the shear direction (Fig. |3]). 



Active liquid crystal with conflicting anchoring 

We now come to the case of an active liquid crystal with the director field along y at y = L 
and along x at y = 0. In the passive case the liquid crystal responds to this conflicting anchor- 
ing by adopting a uniform (in the one elastic constant approximation employed here) splay-bend 
deformation fully contained in the xy plane. 

Fig. |4] shows the director and velocity fields at steady state for small values of the activity 
(. Contrarily to the normal anchoring case, we find no critical non-zero value of for which the 
system is passive. Formally this can be understood as in the passive steady state there are gradients 
of Q which set up a non-trivial shear rate so that the passive solution is always unstable (for C, > 0). 
For small values of C, the active gel expands the region which is aligned along x, while the flow 
field is asymmetric at steady state. 

Fig. Oproves that also in this HAN geometry further 'transitions' are encountered when going 
deeper into the active phase by increasing For > ~ 0.0003, the sign of the flow in steady 
state reverses. Increasing the activity also leads to a more symmetric, Poiseuille-like flow. The 
change in orientation of the fluid flow relative to the director field leads to a strikingly different 
elastic deformation which ends up with a point with ^ = (director completely along the velocity 
gradient direction) within the sample. The director profile can be understood as due to the fluid 
flow in steady state, as it resembles the elastic distortions in steady state found with clockwise 
and anti-clockwise HAN cells under a constant shear flow in Ref. [llSll . It would be of interest to 
develop a semi-analytic theory which accounts for the presence of extra 'transitions', both in the 
HAN and in the normal anchoring geometry, when the activity is increased. 



Dynamic pathway to steady state 

We now consider the dynamic pathway followed by the active liquid crystal to reach its steady 
state. For concreteness, henceforth we will restrict to the case of normal anchoring. 

Fig. |6]show the dynamics followed by a system with L = 40 fxm, which started from a (slightly 
perturbed) spatially constant director field with normal anchoring and no flow. The value of the 
activity parameter is one for which, for L = 10 fim, one ends up with states similar to those in 
Fig. 2. In the larger sample this value of ( is enough to lead to an out-of-plane director field 
in steady state. This is in line with previous observations that larger systems are more easily 
led into an active state [Q]. The early time dynamics is characterised by the generation of a lot 
of deformation in the director field, (reflected in the large spatial changes in Ux^y^z), caused by 
the appearance of spontaneous flow. Thermodynamic forces then set in to reduce the number of 
'bands' or equivalently to reduce the amount of splay and bend. Note also that transiently the 
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director field develops a significant asymmetry, so that it is twisted unevenly, which gradually 
disappears upon approaching the steady state. 

Note that we have also simulated situations in which the activity value is cyclically increased 
and decreased so that the initial condition is different in the two halves of a 'cycle', to compare the 
steady states found in the two cases (data not shown). While with small to moderate values of the 
activity, for which the active phase has just set in, we find no dependence on the initial conditions, 
this is not true as we go deeper in the active phase. In this region we find there is some hysteresis, 
and that the steady state does depend on initial condition for large value of (i.e. deep in the 
active phase). This is perhaps not surprising because when subjected to a flow (in this case stim- 
ulated by the activity), liquid crystals and other viscoelastic gels may exhibit rheochaotic (hence 
strongly initial condition dependent) flow properties [[ItI] . It would be interesting to systematically 
characterise the dependence of steady state and of the dynamic pathway in active gels for large 

and determine whether, under appropriate conditions, complex, or even chaotic, dynamics are 
observed. This is deferred to future work. 



CONCLUSIONS AND DISCUSSIONS 

In conclusions, we have presented a generalised lattice Boltzmann algorithm, inspired by those 
previously used to solve the Beris-Edwards equations of motion of a passive liquid crystal, to in- 
vestigate numerically the hydrodynamic properties of an active (flow-aligning) liquid crystalline 
fluid in an effectively ID geometry. We have studied two distinct boundary conditions (one with 
the director field normally anchored at the boundary plates and another one with conflicting an- 
choring leading to elastic distortions even in the passive phase), and systematically considered 
how the director and flow fields in steady state evolve upon increasing the parameter which is a 
measure of the liquid crystal 'activity'. 

The interplay of activity, which may introduce spontaneous flow in steady state, and the elastic 
distortion caused by the fixed anchoring at the boundaries, gives rise to an interesting phenomenol- 
ogy. We show that our algorithm gives a non-equilibrium transition from a passive to an active 
phase, and that this occurs with zero threshold when there is conflicting anchoring at the bound- 
aries. Deep in the active phase, where semi-analytical methods are typically unhelpful, we find 
further 'transitions' between qualitatively different steady states of flow and director field. When 
normal anchoring is considered, we observe banded flow solutions in which regions of shear- 
oriented filaments or director field are separated by regions of bend-splay distortions (where the 
order parameter decreases slightly), as well as states with director field out of the xz plane in 
steady state. 

The equations of motion we have studied are expected to describe the coarse grained dynamics 
of actomyosin solutions or of microtubular networks with dyneins, when the actin fibers or the 
microtubules have a fixed orientation at the surface (due e.g. to chemical or mechanical treat- 
ment of the surface as in the csae of passive liquid crystal). We note that our algorithm can be 
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straightforwardly applied to 2D and 3D problems; this will be pursued elsewhere. 
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Figure captions 



Fig. 1: Geometry used for the calculations described in the text. The active gel is sandwiched 
between two infinite plates, parallel to the xz plane, lying at y = and y = L. We consider 
(a) normal anchoring and (b) conflicting anchoring (leading to a hybrid aligned nematic cell for a 
passive liquid crystal) 

Fig. 2: Director (left) and velocity (right) fields in steady state for an active liquid crystal 
slab with normal anchoring at the boundaries, and ( equal to 0.0002 (solid line), 0.0004 (long- 
dashed line), 0.0006 (dot-dashed line), and 0.001 (dotted line), respectively. Other parameters are 
as specified in the text. 

Fig. 3: Director (left) and velocity (right) fields in steady state for an active liquid crystal 
slab with normal anchoring at the boundaries, and equal to 0.002 (solid line), 0.003 (long- 
dashed line), 0.005 (dot-dashed line), and 0.007 (dotted line) respectively. Other parameters are 
as specified in the text. 

Fig. 4: Director (left) and velocity (right) fields in steady state for an active liquid crystal 
slab with conflicting anchoring at the boundaries (Fig. [H and equal to 0.00001 (solid line), 
0.00005 (long-dashed line), 0.0001 (dot-dashed line), and 0.0002 (dotted line), respectively. Other 
parameters are as specified in the text. 

Fig. 5: Director (left) and velocity (right) fields in steady state for an active liquid crystal 
slab with conflicting anchoring at the boundaries (Fig. [U and equal to 0.0004 (solid line), 
0.0006 (long-dashed line), 0.0008 (dot-dashed line), and 0.003 (dotted line), respectively. Other 
parameters are as specified in the text. 

Fig. 6: Dynamic pathway to steady state followed by a system with ( = 0.0008, L = 400, and 
other parameters as in the previous section. Solid, long dashed and dot-dashed lines correspond to 
the z, X and y components of the director field respectively. 
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FIG. 1 : Geometry used for the calculations described in the text. The active gel is sandwiched between two 
infinite plates, parallel to the xz plane, lying at y = and y = L. We consider (a) normal anchoring and 
(b) conflicting anchoring (leading to a hybrid aligned nematic cell for a passive liquid crystal). 
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FIG. 2: Director (left) and velocity (right) fields in steady state for an active liquid crystal slab with 
normal anchoring at the boundaries, and C, equal to 0.0002 (solid line), 0.0004 (long-dashed line), 0.0006 
(dot-dashed line), and 0.001 (dotted line), respectively. Other parameters are as specified in the text. 
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FIG. 3: Director (left) and velocity (right) fields in steady state for an active liquid crystal slab with normal 
anchoring at the boundaries, and Q equal to 0.002 (solid line), 0.003 (long-dashed line), 0.005 (dot-dashed 
line), and 0.007 (dotted line) respectively. Other parameters ai^e as specified in the text. 



14 




FIG. 4: Director (left) and velocity (right) fields in steady state for an active liquid crystal slab with con- 
flicting anchoring at the boundaries (Fig. [T] and C equal to 0.00001 (solid line), 0.00005 (long-dashed line), 
0.0001 (dot-dashed line), and 0.0002 (dotted line), respectively. Other parameters are as specified in the text. 
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FIG. 5: Director (left) and velocity (right) fields in steady state for an active liquid crystal slab with 
conflicting anchoring at the boundaries (Fig. [T] and ( equal to 0.0004 (solid hne), 0.0006 (long-dashed 
line), 0.0008 (dot-dashed line), and 0.003 (dotted line), respectively. Other parameters are as specified in 
the text. 
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FIG. 6: Dynamic pathway to steady state followed by a system with ( = 0.0008, L = 400, and other 
parameters as in the previous section.Solid, long dashed and dot-dashed lines correspond to the z, x and y 
components of the director field respectively. 
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